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Abstract. We reanalyze the XMM-Newton data of GRB 011211 showing that the spectral features, interpreted by Reeves et 
al. (2002, 2003) as due thermal emission from a collisionally ionized plasma, can be also reproduced by a reflection model 
(with ionization parameter £ ~ 10 2 ). We discuss the implications of this interpretation, estimating the total mass required in the 
simplified case of a funnel geometry. We conclude that a moderate clumping of the reprocessing material (corresponding to a 
filling factor of the order of / ~ 10~ 3 ) is required. Finally we show that, if this interpretation is correct, a bright quasi-thermal 
component is expected in the optical-UV band (containing about 90% of the luminosity of the illuminating continuum), whose 
presence can be used to test the reflection model. 
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1. Introduction 

The study of spectral features (absorption, emission lines) 
in the X-ray spectra of afterglows of Gamma-Ray Bursts is 
thought to be one of the most powerful tools to probe the cir- 
cumburst environment. In fact, even if the direct observation 
of the central engine is not possible, it might be possible to 
get indirect indication of the composition and structure of the 
GRB environment using emission lines (for recent reviews see 
Lazzati 2002a and Bottcher 2002). Until now, emission features 
have been detected (in some cases with only marginal signifi- 
cance) in 7 GRBs, namely GRB 970508 (Piro et al. 1999), GRB 
970828 (Yoshida et al. 2000), GRB 991216 (Piro et al. 2000), 
GRB 000214 (Antonelli et al. 2000), GRB 011211 (Reeves et 
al. 2002, 2003), GRB 020813 (Butler et al. 2003) and GRB 
030227 (Watson et al. 2003). In most of these cases the X- 
ray spectrum shows the evidence of the Iron line, largely inter- 
preted as due to recombination, either in a photoionized (e.g. 
Vietri et al. 2001) or in a collisionally-ionized (thermal) plasma 
(e.g. Paerels et al. 2000). 

Recently, Reeves et al. (2002, 2003) used a thermal model 
to fit the X-ray spectrum of the afterglow of GRB 011211 ob- 
served by XMM-Newton (unfortunately the observation is af- 
fected by pointing problems and the real significance of the 
lines has been criticized on different grounds by Borozdin & 
Trudolyubov 2002 and Rutledge & Sako 2003). The residuals 
of a fit to the data with a pure absorbed power-law model show 
evident features located at 0.7, 0.9 and 1.2 keV. Reeves et al. 
show that a good fit to the data can be obtained assuming a 
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thermal plasma model (kT ~ 4 keV) enriched in light metals 
(Ca, Mg, Si, S, Ar). The viability of such a model is critically 
discussed by Lazzati (2003) which concludes that a thermal 
origin of the emission puts severe constraints on the parame- 
ters of the emission region: in particular the emitting gas must 
be strongly clumped, with volume filling factors of the order 
of ~ 10~ 6 — 10~ 7 . Another problem comes from the fact that, 
in order to fit the spectrum with the thermal VMEKAL model, 
Reeves et al. are forced to assume overabundance of Mg, S, Si, 
Ar, Ca (about a factor of 10 larger than the solar abundance), 
but a solar abundance of iron, clearly not easily explainable in 
a standard hypernova scenario (Woosley 1993), since the ex- 
ploding star is thought to produce a large amount of Ni, Co and 
Fe (e.g. Woosley & Weaver 1995). 

Reeves et al. (2003) also considered the possibility that the 
spectrum is produced through the reflection by an ionized slab 
(Ross & Fabian 1993; Ballantyne & Ramirez-Ruiz 2001), but 
concluded that this model (either "pure" or with an underly- 
ing power-law) cannot reproduce the data. However (as explic- 
itly noted by the authors), most of the discrepancy between the 
model and the data is localized in the region around 1 keV, 
where lines from elements not included in the ionized reflec- 
tion model, such as S, Ar and Ca, are expected. 

More recently, Butler et al. (2003) and Watson et al. (2003) 
reported the possible detection of emission lines of light el- 
ements (but not Iron) in the Chandra HETGS spectrum of 
GRB 020813 and in the XMM-Newton spectrum of the after- 
glow of GRB 030227. These features can be identified as lines 
of hydrogen-like SXVI and SiXIV (GRB 020813) and from 
helium/hydrogen-like Mg, Si, S, Ar and Ca (GRB 030227). 
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Also in these cases a thermal model could fit the data, but the 
authors noted that these lines can also originate from reflection 
by material illuminated by the burst/afterglow (even though 
specific fits for this case are not presented for GRB 020813). 
Even more interestingly, in the case of GRB 030227 the lines 
are not detected in the first part of the observation (started about 
1 1 h after the burst) but appear only in the last 10 ksec, approx- 
imately 20 h after the trigger. 

All these lines are detected at the ~ 3<r level: although we 
are waiting for a more secure detection, casting away any doubt 
about their reality, the information carried by their presence 
are so important to justify their study. We therefore assume 
that these lines are real, and in this paper we reanalyze the 
XMM-Newton data of GRB 01 121 1 to show that it is possible 
to reproduce them modelling the spectrum with the reflection 
model, once lines from light elements are properly considered. 
We discuss the case in which the reflector has a funnel geome- 
try and we propose simple observational tests for this scenario 
that could be performed with future missions such as SWIFT. 
We find a consistent scenario able to account for the observa- 
tions of the emission lines of light metals (but no iron) in GRB 
011211, and we point out that if this scenario is correct there 
should be an intense optical-UV-soft X-ray emission accom- 
panying the emergence of the emission lines. This is therefore 
a strong test for the reflection interpretation of the X-ray emis- 
sion lines. 

The plane of the paper is the following: in Sect. 2 we report 
the results of the fit with the reflection model. In Sect. 3 we 
derive and discuss the constraints on the GRB environment and 
the properties of the reflection model. 

In the paper the following cosmology is assumed: Hq = 70 
km s- 1 Mpc- 1 n A = 0.7 and Q=l. 

2. Spectral fits 

XMM-Newton data of GRB 011211 processed through the 
standard pipeline are publicly available We extract the PN 
spectrum of the first 5 ks, where lines are visible (Reeves et 
al. 2002, 2003), using the standard package SAS. We fitted the 
spectrum using XSPEC (v.ll). 

2. 1. Pure reflection 

Unfortunately, the reflection model available in XSPEC (avail- 
able as table model file; see the description of the model in 
Ballantyne, Iwasawa & Fabian 2001) is calculated for two fixed 
values of the abundance, solar and 2 x solar. Moreover, as antic- 
ipated in the Introduction, it does not include light metals such 
as S, Ar and Ca. The reflection model assumes that a slab is 
illuminated by a power law continuum and the basic fitting pa- 
rameter is the ionization parameter, defined as £ = L\\\j (nR 2 ), 
where Lu\ is the luminosity of the continuum, n is the density 
of the illuminated slab and R is the distance between the slab 
and the illuminating source. 

To include the contribution of light metals we fit the data 
using the following procedure: we start fitting the data by us- 

1 http://xmm.vilspa.esa.es/external/xrnm_sched/too 
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Fig. 1. Spectrum (top panel) and data/model ratio (bottom 
panel) of the first 5 ksec of the PN data fitted with the reflec- 
tion model and lines from light elements (see discussion in the 
text). 



ing the (absorbed with Galactic iV H = 4.2 x 10 ' cm , 
Reeves et al. 2003) reflection model with free redshift. In the 
fitting procedure we excluded all the bins below 0.3 keV, due 
to known problems with the calibration of the response ma- 
trix (e.g. Brinkmann et al. 2001). As expected, the residuals 
of the best fit (which converges at z = 1.9) show clear devi- 
ations around 0.7, 0.9, 1.1, 1.2 and 1.5 keV. This agrees with 
the analysis by Reeves et al. The line at 0.7 keV can be at- 
tributed to SiXIV. The model includes this element, but the fit 
of the line is poor, since the data requires a brighter line (which 
may correspond to a larger abundance). The other features can 
be identified as emission lines from elements not included in 
the reflection model (SXVI at 0.9 keV, ArXVII and SXV at 
1.1 keV, Ar XVIII at 1.2 keV and CaXX at 1.5 keV). We then 
added Gaussian emission lines to account for these residuals. 
The folded spectrum and the model/data ratio are shown in 
Fig. 1. The unfolded spectrum is shown on Fig. 2. The best fit 
(x 2 /d.o.f. = 27.45/26) parameters are: Logi = 2.25 ± 0.30, 
z = 1.93 ±0.10 with an incident power-law of photon spectral 
index T = 2.43 ± 0.20. 

The parameters of the lines we have added, but that are not 
included in the fitting model, are reported in Table 1 . Similarly 
to Reeves et al. (2003) we find that the lines require a red- 
shift z = 1.8-1.9, lower than the redshift of the GRB, 
z = 2.140 ± 0.001 (Holland et al. 2002), suggesting that the 
emitting material is moving with a velocity v ~ 0.1c. The lu- 
minosity of the Fe line, included in the model is 2.3 x 10 44 
erg s — 1 . For the possible presence of the H-like Nickel (rest 
frame energy E = 8.102 keV) we obtained an upper limit to 
the equivalent width of EW < 450 eV (Reeves et al. 2003 
found EW < 800 eV). The expected luminosities of light el- 
ements, compared to iron, have been calculated by Lazzati et 
al. (2002b). The maximum contrast between the luminosity of 
the lines of the light metals and the iron line is reached around 
£ ~ 10 2 . For £ ~ 2 x 10 2 , required by our fit, we expect a 
luminosity ratio of the order of 3, in agreement with the lumi- 
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Fig. 2. Unfolded spectrum of the PN spectrum fitted with the 
reflection model. The narrow lines at 0.7, 0.9, 1.1, 1.2, and 1,5 
keV are the lines of the elements reported in Table 1, not in- 
cluded in the reflection model. 
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Table 1. Parameters of the lines added to the reflection model 
(see text for details). Luminosities refer to the pure reflection 
case; in the reflection+power-law case they have to scaled by a 
factor 0.7. 



nosities reported in Table 1 . Looking in detail at each element 
we note that lines of S and Ar appear to be more luminous than 
expected, by a factor of 2 and 5, respectively (note however 
that, due to the relatively large uncertainties, the normalizations 
are consistent with solar abundance within ler and 2a, respec- 
tively). This discrepancy (present also in the fits presented by 
Reeves at al.) could be ascribed to a larger abundance of these 
elements with respect to the solar value. 

Thus we can conclude that, once light metals are included, 
the reflection model provides a consistent fit of the XMM- 
Newton spectrum. 

2.2. Reflection plus underlying power-law 

Although a pure reflection model can account for the observed 
X-ray spectrum, it is possible, and even likely, that part of the 
flux originates from the afterglow continuum. It is important 
to stress that, unlike the case of AGNs, the continuum repro- 
cessed through the reflection and that directly observed from 
the afterglow, simultaneously with the lines, could be differ- 
ent. In fact the radiation reprocessed by the reflector originates 



from the unobserved, rapidly changing burst/early-afterglow 
phase (hereafter illuminating continuum), while the afterglow 
radiation (observed when lines are detected) is produced much 
later (hereafter late afterglow continuum). 

We performed a series of fits to asses how much of the total 
emission can be ascribed to an underlying power-law contin- 
uum produced by the late afterglow continuum . We fit the data 
with a model including a power-law and the reflection spec- 
trum. We fix the relative normalization of the reflection spec- 
trum and of the gaussian lines to the value found in the fit with 
the pure reflection case. 

First we constrain the two slopes, that of the illuminating 
continuum and that of the late afterglow power-law, to the same 
value. The fit converges to T = 2.26 ± 0.1, Log£ = 1.9 ± 
0.1 with x 2 /d-o.f. = 23.71/31. The significantly improves 
with respect to the pure reflection case. The presence of a direct 
power-law is therefore suggested by the data. 

We then allowed the two photon indices to vary inde- 
pendently. In this case the best fit parameters are: T p i = 
2.04+°;^, r refl = 2.39l°;J 6 , Logti = 2.1 ± 0.1 with 
X 2 /d.o.f. = 23.13/30. Within 90% conf. level, the two slopes 
are consistent. Reflection contributes for ~70% of the flux. 
Correspondingly, the luminosity of the emission lines given in 
Tab. 1 decreses by a factor 0.7. 

3. Discussion 

We showed that the XMM-Newton PN data of GRB 011211 
can be satisfactorily fitted by the pure reflection model, if lines 
from the elements not included in the spectral model are prop- 
erly considered. The ionization parameter £ is well constrained 
to be of the order of 10 2 , since the data do not present a promi- 
nent iron line while lines of light metals are luminous. As men- 
tioned, the absence of a prominent iron line is naturally pre- 
dicted when the ionization parameter £ is around 10 2 (Lazzati 
et al. 2002b). In this case, due to the importance of the Auger 
effect, the emission from iron is strongly depressed compared 
to the other metals. 

The reflection model can explain very different spectra, 
since variations of £ in different bursts could produce different 
reflection spectra, with prominent or almost absent lines (e.g. 
Ballantyne & Ramirez-Ruiz 2001). Low values of the ioniza- 
tion parameter (£ < 100) will produce a reflection component 
with evident lines. Large ionization parameter (£ > 10 4 ), on 
the other hand, would produce almost featureless spectra, with 
a shape similar to the shape of the illuminating continuum. For 
a narrow range of £ centered around 10 2 the spectra will show 
luminous lines from light metals and a depressed K Q iron line. 
It may be argued that in the case of GRBs a well defined ion- 
ization parameter £ can hardly be attained. In the geometry we 
consider here the observer receives a spectrum resulting from 
the convolution of reflected components with different illumi- 
nating spectrum. Regions closer to the line of sight will con- 
tribute low £ components, while regions away from the line of 
sight will contribute higher £ components, since the time-delay 
is longer and the ionizing continuum is a monotonically de- 
creasing function (at least in the afterglow phase). While this 
is true, it is clear that one of the components (and hence a well 
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defined £) will dominate. Let us now consider the early X-ray 
afterglow as illuminator, since the GRB proper will have a too 
high £ to produce lines, as discussed above. If the illumina- 
tor has a power-law decay in time L oc t~ a , the more (less) 
luminous phase will dominate the observed reflected compo- 
nent for a > 1 (a < 1). Unfortunately the X-ray afterglow 
has never been observed early enough to constrain its decay 
in the early phases, but it is not unreasonable to assume that 
the X-ray band lies initially below the cooling frequency, and 
that its decay slope is initially the same as the optical slope at 
later stages. In this case, bursts with an initially shallow decay 
in the optical, such as GRB 011211 (Holland et al. 2002) and 
GRB 020813 (Covino et al. 2003), should show low-ionization 
features, while bursts with initial fast decay such GRB 991216 
(Halpern et al. 2000) should have X-ray spectra with strong 
iron lines. Even though it is not possible to make statistics with 
such a handful of cases, it is tantalizing that such a correlation 
is indeed respected in all cases. 

In the following we derive the geometrical and physical 
setup of the reflecting material. We note that the actual geom- 
etry of the funnel is unimportant for most of the estimates re- 
ported below, aimed to find some basic quantities such as the 
overall size and the total mass of the line emitting material. 
The emitting layer is supposed to be illuminated by an intense 
X-ray flux: this requires that either the illuminating source is 
less collimated than the funnel (if it has a conical profile), or 
that the funnel deviates from a perfect conical geometry, being 
e.g. parabolic. In these conditions the funnel walls can receive 
a non-negligible fraction of the ionizing flux. We then consider 
more complex geometries, involving clumps of overdense ma- 
terial, possibly located close to the funnel surface, as illustrated 
in Fig.|3] 

3.1. General constraints 

Once the ionization parameter and the line luminosities are 
known, we can infer important clues on the physical state of 
the reflector. The "observed" or fixed quantities are the follow- 
ing: 

1. the luminosity of the line L\[ ne ; 

2. its duration tii ne ; 

3. the time i app at which the line starts to be visible with the 
luminosity L\i UC ; 

4. the ionization parameter £; 

5. the efficiency r\ of converting the illuminating continuum 
into the specific line; 

6. the total mass Mr (which we require not to exceed 
~1OM ). 

Note that the time at which lines appear (t app ) is an observable 
in principle, but it may require early observations of the X-ray 
afterglow to be detected. Indeed, it has been already determined 
only in one case, (GRB 030227) while we have only an upper 
limit for all other cases. However, it is likely that with a short 
pointing time, like the one provided by SWIFT, the number of 
cases in which this parameter will be observable will strongly 
increase. Therefore, in view of a general discussion, we sup- 
pose that t app is known. For the specific case of GRB 011211 




Fig. 3. Sketch of the geometry assumed for the reflecting ma- 
terial (not to scale). The material is organized in clumps along 
the walls of the funnel (assumed to have a parabolic shape). 
Clumps could also lie along the line of sight. 

we only have an upper limit to t app , since the lines were de- 
tected in the first part of the observation. 

Unknown physical quantities of the system are: 

1 . the distance R between the illuminator and the line emitting 
material. 

2. the collimation angle 9 of the illuminator (see Fig. 3). 

3. the electron density n c (which, for the large ionization con- 
ditions discussed here is close to the proton density n p ); 

4. the illuminating energy Em (or, alternatively, the illuminat- 
ing luminosity Liu) directly intercepted by the line-emitting 
material. Part of this illuminating energy can come from 
the burst emission, and part from the early afterglow. It is a 
rapidly variable quantity, and for the sake of simplicity we 
will assume that Lm is an appropriate average over tin; 

5. the time im for which the illumination lasts; 

6. the number of lines photons, k, produced by each atom dur- 
ing tin; 

7. the filling factor /; 

8. the height of the shell Ai?; 

Thus we have a total of 6 quantities provided by observa- 
tions and 8 unknown quantities. Therefore we will consider 9 
and ti\\ as free parameters, even if their values are in any case 
limited (for instance, "reasonable" aperture angles of the illu- 
minator should be between 10 and 60 degrees; the illuminating 
time should be shorter than the time for which the line is visi- 
ble, till < tline)- 
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In the following we outline how we can connect these two 
sets of parameters and how we can derive the parameters of the 
system by the observational quantities. All the time intervals 
(hlh tu n e, t app ) are intended as measured in the frame of the 
burst. We adopt the notation Q = 10 X Q X and use cgs units. 

The size of the reflector is related to the time at which lines 
become visible: 



R = 



ct 



app 



1 — COS ( 



(1) 



The total mass of the element emitting a specific line (denoted 
by the subscript i) is related to the number TVj of photons in 
that line by: 



NiWim v 



(2) 



where Wi is the atomic weight of the element, k is the number 
of line photons a single atom can emit during the illumination 
phase (Lazzati et al. 1999): 



k = 



tree 



where the recombination time is: 
1 

^rcc — 8 

a v n c 

n e is the numerical density of the free electrons and a r is: 

0.4961 



5.2 x W-^ZX 1 / 2 



0.429 + 0.5 ln(A) 



AV3 



(3) 



(4) 



(5) 



where A = 1.58 x 10 5 Z 2 T~ 1 , Z the atomic number of the 
element and T the electron temperature. 

The total number of photons produced by atoms of the ele- 
ment i can be written as: 



Ni = 



-^linc^linc 



(6) 



where En nc = Limc^imc is the total energy in the line and 
is the energy of the line. The mass of the layer emitting the 
observed line is (using the mass abundance A.- L of the emitting 
element): 



M ' " A = 



Mi _ Ln nc tn ne W l m p 



(7) 



This mass is contained into a layer of volume V, of thickness 
Ax, radius R sin 9 and height AR. The material reflecting the 
illuminating continuum and producing the line is contained in 
a layer of height Air of unit scattering optical depth: tt = 
<7Tn e Ax = 1. Since the matter is highly ionized the density 
n e of free electrons is the same as the proton density n p . This 
allows us to derive the mass of this layer in a way independent 
from Eq. @ and also independent of the electron density: 



Mr = 



n p m p V = 2nRsm9ARAxn e m p 



2tt 

— RAR sin 6m p 
57ri?Ai?sin6> g 



(8) 



We then have two independent estimates of the mass involved 
in the line production, Eq. and Eq. (jSJi. By equating them 
we can derive the value of k (which becomes independent of 



_ W l m p 



Aiei 5wRARsint 



(9) 



From the definition (|3} and Eq. @ and inserting k from l|9} we 
can derive the density: 



k 



a r tm 

WiTUp 



-^linc^lii 



cm 



AiCi 5nRAR sin 9a T tm 
The ionization parameter is defined as: 



n c R 2 



Inserting the density given by Eq. dlOi . we have: 



1 R Wim p E hnc _ x 
■ £ -r-^ — r—^ : — erg s 



5tt AR Ai€i sin(9a r iin 
The efficiency 77 in producing the line is defined as: 

-^linc -^linc^line 

Em iilliill 

Using Eq. il 2i and Eq. Jl 3i we derive: 

AR _ 1 WiTUp 
R 5n ^ AiEi sin 9a r 

From Eq. © the mass of the layer turns out to be: 

Mx = R\r l Wim P - 



Ai6ia r 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



The total amount of mass M* (i.e. contained in a shell of sur- 
face AttR 2 and height AR and assuming the same density) will 
be: 



,AR 



(16) 



We impose that this mass cannot exceed Mt ~ 10M Q . If M* 
exceeds this value, we assume that the material is distributed 
inhomogeneously, more dense in the walls of the funnel or or- 
ganized in single clumps (see below). We parametrize the de- 
viation from the homogenoeus case by the filling factor: 



f = Ml 



(17) 



Note that there is in principle another equation relating AR 
with the duration of the line tu ne given by 



AR 



ctm 



cos 9 



(18) 



This equation is valid only in the case in which the duration 
of the emission of the line is dominated by geometrical effects, 
i.e. tin < AR/c. In principle, in this case it is possible to fix 
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the aperture angle of the funnel, using Eq. dl4> together with 
Eq. flBJ and Eq. (Q): 

tan0(l-cos0) ~ — ^rjEiH^h^ (l9) 

57T AiCia T ilino 

For the application to the specific case of GRB 01 121 1 we have 
preferred to use the angle 6 as a free parameter, using Eq. dl 81 
as a consistency check. We do that because of the uncertainties 
related to the specific geometry. 

We stress that the above estimates have been derived for a 
specific geometry, i.e. a uniform funnel over-dense with respect 
to the rest of the remnant. It is conceivable that the real situation 
is much more complex: for instance, dense clumps of matter 
could be embedded in a less dense medium, and be located also 
in the interior of the funnel, not only on its walls (Fig.0. On 
the other hand the values of the parameters we derive will not 
dramatically change in this case, and can therefore be thought 
as indicative numbers. 

3.2. Application to GRB 011211 

As an illustrative example of the application of the estimates 
given above we can apply our treatment to the specific case of 
GRB 011211. For the observational quantities we set: 

- L Une = 7 x 10 44 erg s- 1 for the SiXIV line; 

- £iinc = 5x 10 3 /(1 + z) ~ 1.7 x 10 3 seconds; 

- t app < 4 x 10 4 /(1 + z) ~ 1.3 x 10 4 seconds; 

- £ = io 2 6; 

- rj = 10~ 2 ?7_2 (see Ghisellini et al. 2002); 

- M T = 10M Q 

In addition, we have used the redshift z = 2.14; the line en- 
ergy esixiv — 1-99 keV and an abundance A; equal to twice 
the solar value. The temperature of a photoionized plasma illu- 
minated with £ ~ 100 is predicted to be in the range 10 5 — 10 6 
K (e.g. Kallman & McCray 1982). We have then assumed that 
the line emitting material is at T = 5 x 10 5 K. A variation of 
the temperature, directly affecting the recombination rate, im- 
plies a different amount of matter necessary to account for the 
emission and therefore the estimate of the density and of the 
clumping factor. Eq.Q indicates that the dependence of t Iec 
on the temperature is t Iec ~ T 3 / 4 . A larger (lower) tempera- 
ture would decrease (increase) the recombination rate, resulting 
in a larger (lower) clumping / ~ T~ 3 / 4 . 

As discussed in Sect. 3.1, in the specific case of GRB 
011211, since the appearence time of the line is unknown, the 
value of t app given above should be considered an upper limit 
to the actual value and, therefore, the estimates given below 
should be considered as limits. In particular the derived dis- 
tance R is un upper limit, the density is a lower limit and thus 
the clumping factor should be considered as a lower limit. In 
Fig. 4 we show a few derived quantities as a function of the 
aperture angle 9 and for txxx — 300 and 1000 seconds. The dis- 
tance R is of the order of a few times 10 15 cm; the density n e 
is of the order of 10 14 -10 15 cm~ 3 ; AR/R is of the order of a 
tenth, and the filling factor of the order of / ~ 10~ 3 . The mass 
of the line emitting layer is M\ = 10~ 4 -10~ 2 Mq. It is inter- 
esting to note that the assumed (admittedly simple) geometry 



-2 



1 1 


1 1 1 

t u 


1 1 1 1 1 1 1 1 1 1 

,= 103 s 


i i i 1 i i i i 














AR/R 






1 1 


1 1 


i i i i i i i i 1 i 


1 1 1 1 IJ— 






1= 300 s 


n u : 










1 1 


1 1 1 


i i i 1 i i i i 1 i 


1 1 1 1 1 1 1 1 



10 20 30 40 50 60 
6 [degrees] 

Fig. 4. Derived quantities for GRB 02121 1, for t m = 10 3 s (top 
panel) and for tux — 300 s (bottom panel). The filling factor is 
in units of 10 -3 , the distance R in units of 10 15 cm, the density 
n in units of 10 14 cm -3 and the mass of the emitting layer in 
units of 10 -3 solar masses. The quantities R, AR/R and M\ 
do not depend on tax. The dashed line in the top panel is the 
value of AR/R derived using Eq. Q and Eq. Jl 81 . See text for 
the assumed values of the input parameters. 

can also account for the line duration, which yield an indepen- 
dent value of AR/R (dashed line in the top panel of Fig. 4). 
According to it, the illumination angle turns out to be ~ 45° 
(i.e. where the dashed and solid lines labeled as AR/R cross). 

The portion Em of the total energy of the illuminating con- 
tinuum directly impinging on the reflecting matter will be 

E m ~ -L linc t linc 3.5 x 10 50 erg. (20) 

The total apparent isotropic energy contained in the illumi- 
nating component will be Ei SO = Exxi/ACl, where Ail is the 
solid angle subtended by the reflecting matter as observed the 
location of the illuminating source. It can be shown that in the 
geometry assumed here and with the geometrical parameters 
found above Afl ~ 0.01 4n, giving E{ so ~ 3 x 10 51 erg, con- 
sistent with the early phase of the afterglow. 

We conclude that the reflection model offers a good inter- 
pretation of the data. The reflecting material must be clumped, 
although the required clumping factor is smaller (and the cor- 
responding filling factor greater) than the one required by the 
thermal model (Lazzati 2003). 

3.3. Re-emission of absorbed radiation 

An important consequence of the assumption that reflection is 
responsible for the production of detected lines is that, espe- 
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cially when £ is low (£ < 10 2 ) as in the case discussed here, 
most of the incoming radiation is absorbed by the reflecting 
material and successively re-emitted as low-frequency radia- 
tion. 

In the low-ionization condition discussed here the reflected 
luminosity is ~ 10% of the total incident luminosity (e.g. Zicky 
et al. 1994): therefore the remaining ~ 90% is absorbed by 
the medium, which is heated-up and re-emits what it has ab- 
sorbed. The luminosity L rcpl of this reprocessed component is 
rather important and amounts to: 

t n n -^line n i r\46 -^linc,45 — 1 /m\ 

£repr = 0.9 ~ 9 x 10 ■ — erg s (21) 

V V-2 

This emission lasts exactly for the same time as the X-ray 
emission line, and is therefore a powerful diagnostic of the 
reflection model. The above luminosity, for the case of GRB 
011211 at z = 2.14, corresponds to a flux of 3 x 10~ 12 erg s _1 
cm -2 . The problem is to assess the spectral region where this 
emission peaks. We should solve the radiation transfer into a 
slab, heated by the incoming rapidly variable high energy emis- 
sion from the early afterglow. The surface of the slab, directly 
illuminated by the incoming continuum will be heated to large 
temperatures, close to the Compton temperature Tc ~ 10 7 K. 
The emission from these outer portion will peak in the X-ray 
region, at energies around 1 keV. Gas in the layers deeper into 
the slab will be close to the thermal equilibrium: therefore the 
temperature of these region will be close to the black-body 
temperature corresponding to a black-body emitting a lumi- 
nosity of the order of Lbb = ^repr^ 

/ r \ 1/4 / t \ 1/4 

where the area A = 2nRs\\\6/S.R has been assumed of the 
order of A — 10 30 A3o cm 2 . The maximum of the emission 
falls at the frequency v hh = 3fcT/[(l + z)h] ~ 10 15 Hz, i.e. in 
the near UV. 

The actual outcoming spectrum will be a complex integral 
over the emission from the different layers (weighted by the 
optical depth of the layer), in strongly time-dependent condi- 
tions. The detailed calculation of the spectrum is beyond the 
scope of this paper. Moreover, note that the conditions appro- 
priate to analyze even a steady state case are at the limit of 
the calculation possibility of the program CLOUDY (Ferland 
1996) which gives results which are not completely reliable. 

What we can safely conclude is that, simultaneously with 
the presence of the X-ray emission lines from light metals, 
there should be a bright component between the UV and the 
soft X-rays, which has a luminosity ^100 times larger that the 
line luminosities (and 10 times larger than the entire X-ray re- 
flection continuum). The low frequency tail of this reprocessed 
component falls in the optical band, where can be observed as a 
rebrightening of the light curve. As an illustrative example, if 1 
per cent of the reprocessed luminosity is emitted in the optical 
V-band, we should see a magnitude of V ~ 22. 

We stress that the presence of the reprocessed emission is 
a natural consequence of the reflection interpretation: obser- 
vations in the UV-soft X-ray bands, confirming or ruling out 



the existence of this component accompanying the X-ray emis- 
sion lines, provide a powerful and reliable diagnostic tool to 
study the origin of the emission features. The only way we en- 
visage not to see this reprocessed component in the reflection 
scenario (with a small ionization parameter) is to have dense 
and small clumps with a Thomson optical depth close to unity. 
These blobs would efficiently reflect the incoming continuum, 
but only a small fraction of it is absorbed and re-emitted. Since 
this requires a rather ad-hoc setup, we consider it unlikely. 

3.4. Close reprocessor 

We discussed the case of a far reprocessor, located at a distance 
R ~ 10 15 — 10 16 cm from the illuminator. There is however 
the possibility that the material producing the lines is located 
very close to the central source. In this case the duration of the 
line is linked to the time for wich the central engine is active. 
A model of this type has been proposed by Rees & Meszaros 
(2000) and Meszaros & Rees (2001). 

The discussion given in the previous Sections can be eas- 
ily applied to this different scenario. In this case Eq. Q 

can 

be dropped, since the geometrical effects are not important in 
determining the time at which the lines appear. All the other 
equations still apply, once 9 is fixed to tt/2. In this way we can 
model a cilindrical geometry, with the illuminating source lo- 
cated on the axis. We fixed the size to R = 10 13 cm, the typical 
size considered by Rees & Meszaros (2000) and Meszaros & 
Rees (2001). 

We can derive the interesting physical parameters for the 
value of the observational quantities given above in the same 
way as above. We found: / = 4, n e = 1.2 x 10 20 cm~ 3 , 
AR/R = 0.03 and M t = 2.5 x 10~ 8 Mq. The large value of 
the density is a consequence of the fact that, for a fixed value 
of the ionization parameter, Eq. (II Tj implies a large density for 
compact (i.e. small R) reprocessor. The value of / (close to 
1) shows that no clumping is needed in this case, due to the 
smaller quantity of iron required compared to the far reproces- 
sor case. 

Although the values of the derived physical quantities are 
reasonable, this version of the reflection scenario suffers of 
many problems. First of all (e.g. Lazzati et al. 2002b), the lines 
and the continuum responsible for that lines are observed to- 
gether: thus it is not possible to produce lines with the large 
observed EW (hundreds of eV). This problem is naturally 
by-passed by the far-reprocessor scenario, since, due to the 
crossing time delay, the lines are observed much later than the 
direct continuum. A possibility to solve the difficulty in the 
close reprocessor is to admit that the illuminating continuum is 
strongly anisotropic and threfore it is not directly observable. 

Another problem is related to the large temperature reached 
by the illuminated material. In fact, following the results of 
Sect. 3.3, the temperature of the material will lie between Tbb 
and Tc- Due to the small area, A ~ 10 25 cm 2 , Eq. (1221 gives 
Thb ~ 3.5 x 10 6 K. The reprocessed thermal radiation will 
therefore peak in the soft X-ray domain (> 0.1 keV). A strong 
soft excess has never been observed in the X-ray spectra of af- 
terglows with lines. 
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For these reasons we consider this possibility unlikely. 

3.5. Time evolution 

All the issues related to the GRB are characterized by a strong 
time-dependent phenomenology. Besides the normal afterglow 
behavior produced by the decelerating fireball, in the X-ray 
band we have the possibility to see the reflection component, 
which it will too be variable in flux and possibly in shape. The 
visibility of the reflected component with the associated emis- 
sion lines depends upon the relative strength of the afterglow 
and of the reflected component, and on the ionization param- 
eter determining which lines are most efficiently produced. In 
general, the reflection component will be visible (if at all) only 
in intermediate phases of the afterglow: on one hand the lumi- 
nosity of the afterglow must decrease sufficiently to allow the 
reflection to be observed, but on the other hand the reflected 
radiation can be observed only for a time smaller than the light 
travel time of the line emitting material (after that time we ex- 
pect the reflection component to disappear). 

At these intermediate times, if the reflection component 
contributes to the total X-ray flux, it will produce a "bump" 
in the X-ray lightcurve, flagging the likely presence of lines 
in the X-ray spectrum. At the same time, we expect a change 
in the slope of the X-ray spectrum, which will depend upon 
the ionization parameter. In fact, for relatively small £, the re- 
flecting layer will not be completely ionized, and the reflection 
component will be hard (i.e. characterized by a photon spec- 
tral index T < 2) in the 0.1-10 keV band, while for large £ 
the reflecting layer will be almost completely ionized, and the 
shape of the reflected component would closely resemble the 
shape of the illuminating continuum. These simple arguments 
offer a direct test of the reflection scenario: suppose that the 
X-ray light curve, after a few hours, flattens: if this flattening 
is due to the emergence of the reflection component then its 
spectrum should be harder or almost equal to the spectrum of 
the "normal" afterglow. If it is harder we expect a small ion- 
ization parameter and the presence of emission lines from light 
metals (and weak or absent lines from iron, cobalt or nickel). 
If the spectrum maintains a constant or quasi-constant shape, 
then we expect the presence of K Q iron (or cobalt o nickel) lines 
and no emission from the fully ionized light metals, unless the 
ionization parameter is so large to maintain also the iron fully 
ionized. 

In conclusion we have presented a coherent scenario in 
which the emission of lines from light metals can be explained 
as due to the reflection by dense matter around the GRB. This 
scenario can be tested by observations. In particular the predic- 
tion of a bright optical-UV component coming together with 
the appearence of the emission lines offers a strong test to the 
overall picture. 
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